Quantum phase transitions in bilayer quantum Hall systems 
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We construct a quantum Ginsburg-Laudau theory to study the quantum phase transition from 
the excitonic superfluid (ESF) to a possible pseudo-spin density wave (PSDW) at some intermediate 
distances driven by the magneto-roton minimum collapsing at a finite wavevector. We analyze the 
properties of the PSDW and explicitly show that a square lattice is the favorite lattice. We suggest 
that correlated hopping of vacancies in the active and passive layers in the PSDW state leads to 
very large and temperature dependent drag consistent with the experimental data. Comparisons 
with previous microscopic numerical calculations are made. Further experimental implications are 
given. 



Introduction. Spin-polarized Bilayer Quantum Hall 
systems at total filling factor vt = ^ have been un- 
der enormous experimental and theoretical investigations 
over the last decade P, Q . When the interlayer separa- 
tion d is sufficiently large, the bilayer system decouples 
into two separate compressible v — 1/2 Fermi Liquid 
(FL) layers. However, when d is smaller than a crit- 
ical distance, even in the absence of interlayer tunnel- 
ing, the system undergoes a quantum phase transition 
into a novel spontaneous interlayer coherent incompress- 
ible phase which is an excitonic superfluid state (ESF) 
in the pseudospin channel 0, [3, 0, Although there 
are very little dissipations in both the ESF and FL, the 
experiment ^ discovered strong enhancement of drag 
and dissipations in an intermediate distance regime. One 
of the outstanding problems in BLQH is to understand 
novel phases and quantum phase transitions as the dis- 
tance between the two layers is changed. If the exper- 
imental observations indicate that there is an interme- 
diate phase between the two phases remains controver- 
sial. Even it exists, different scenarios are proposed for 
the nature of the intermediate phase (1, 0, B Q • Using 
Hartree-Fock (HF) in the Lowest Landau Level (LLL) 
or trial wavefunctions approximations, many authors |6[ 
proposed different kinds of translational symmetry break- 
ing ground states as candidates of the intermediate state. 
In this paper, we develop a quantum Ginsburg-Landau 
(QGL) theory [13] to study the nature of the intermedi- 
ate phase. We propose there are two critical distances 
del < dc2 and three phases as the distance increases. 
When < d < dd, the system is in the phase ordered 
ESF state which breaks the internal U{1) symmetry, 
when del < d < dc2, the system is in a pseudo-spin den- 
sity wave ( PSDW ) state which breaks the translational 
symmetry, there is a first order transition at dd driven by 
magneto-roton minimum collapsing at a finite wavevec- 
tor in the pseudo-spin channel. When c?c2 < d < oo, the 
system becomes two weakly coupled v = 1/2 Composite 
Fermion Fermi Liquid ( FL) state. There is also a first 
order transition at d = dc2- However, disorders could 
smear the two first order transitions into two second order 
transitions. We construct a quantum Ginzburg-Landau 



theory to describe the ESF to the PSDW which break 
the two completely different symmetries and analyze in 
detail the properties of the PSDW. By using the QGL, 
we explicitly show that a square lattice is the favorite 
lattice. The correlated hopping of vacancies in the active 
and passive layers in the PSDW state leads to very large 
and temperature dependent drag consistent with the ex- 
perimental data in [5|. Recently, the effects of small im- 
balance above the PSDW were studied in pT| and found 
to explain the experimental observations nicely in 12l |. 



Formalism in phase representation. Consider a 
bi-layer system with Ni { N2 ) electrons in top ( bottom 
) layer and with interlayer distance d in the presence of 
magnetic field B = V x A: 
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where electrons have hare mass m and carry charge — e, 
Ca,a — 1,2 are electron operators in top and bottom 



layers, Spa{x) 



= 4(x)c„(f) 



— 1,2 are normal 



ordered electron densities on each layer. The intralayer 
interactions are Vn — V22 = e^/er, while interlayer in- 
teraction is V12 = V21 — e^/eVr^ -|- d^ where e is the 
dielectric constant. For simplicity, we only discuss the 
balanced case in this paper. The effects of imbalance 
were discussed in detail in lol. 11 1. 



Performing a singular gauge transformation 4>aix) = 
g^/d=a:'ar^(^-^^')p(^')c^(f) where p(f) = c\ix)ci{x) + 
C2{x)c2{x) is the total density of the bi-layer system. We 
can transform the Hamiltonian Eqn[T]into a Lagrangian 
of the Composite Boson (jja coupled to a Chern-Simon 
gauge field [loj . We can write the two bosons in terms 
of magnitude and phase (f)a = y/pa + Spae^^" , then after 
absorbing the external gauge potential A into a, we get 
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the Lagrangian in the Coulomb gauge : 

L = z5p+(ia.0+-ao)+2^[iv0+~a]2 
1 i 

where i5/9± = Spi±Sp2, 0± = 6i±62, they satisfy commu- 
tation relations [5pa{x), 9i3{x')\ = 2iHSa0S{x — x'), a,(3 — 

±. -p = p^ + p^^v± = y^^. 

It was shown in [lo| that the functional form in the spin 
sector in Eqn [2] achieved from the CB theory is the same 
as the microscopic LLL+HF approximation achieved in 
• However, the spin stiffness and the V- (q) in Eqn[2] 
should be replaced by the effective ones calculated by 
the LLL+HF approximation: pE and Ve((z) = a — bq + 
cq^ where the non-analytic term is due to the long-range 
Coulomb interaction, a ^ (P,b ~ cf, but c remains a 
constant at small distances. In the ESF state it is 
convenient to integrate out Sp- in favor of the phase field 



^s = ^ildr9-r+p4yo_r (3) 

where the dispersion relation of the NGM including 
higher orders of momentum can be extracted: 



= [2pEVE{q)]q'^ :^q^{a-bq + cq^ 



(4) 



Instability driven by the collapsing of magneto- 
roton minimum. As shown in [16| in the context of 
possible supersolid in Helium 4, the advantage to extend 
the dispersion relation beyond the leading order is that 
the QGL action can even capture possible phase transi- 
tions between competing orders due to competing inter- 
actions on microscopic length scales. In the following, we 
study the instability of the ESF state as the distance in- 
creases. It can be shown that the dispersion curve Eqn[4] 
and VE{q) indeed has the shape shown in Fig. la [l^ and 
Fig. lb respectively. 

The phenomenon of the collapsing of the magento- 
roton minimum as the distance increases was clearly 
demonstrated in the numerical calculations in the 
LLL4-HF approaches [1, 14 1 and detected by inelastic 
light scattering [l^. It was estimated that q^l ~ 1, so 
the lattice constant of the PSDW is of the same order of 
magnetic length I which is ~ lOOA. The critical distance 
del is also of the same order of the magnetic length. In re- 
ality, the instability happens before the minimum touches 
zero. 

Effective action in the dual density representa- 
tion and a Feymann relation in the pseudo-spin 
channel. Because the original instability comes from the 
density-density interaction, it is convenient to integrate 
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FIG. 1: The zero temperature phase diagram in the balanced 
case as the distance between the two layers increases. ESF 
where < i/> >^ 0, < >= stands for excitonic superfluid, 
PSDW where < ip >= 0, < stands for pseudo- 

spin density wave phase, FL stands for Fermi Liquid, (a) 
Energy dispersion relation 01(9) in these phases, (b) VE{q) 
in these phases. The cross in the PSDW means the negative 
minimum value of Ve (q) is replaced by the PSDW. The order 
parameters are also shown. In fact, the instability happens 
before the minimum touches zero. 



out the phase field in favor of the density operator. Ne- 
glecting the vortex excitations in 0_ and integrating out 
the 9- in EqnO leads to: 



C[Sp^] = -5p-{- 



2pEq^ 



VE{q)]Sp iq,ujn) 



(5) 

where we can identify the dynamic pseudo-spin 
density-density correlation function S-{q,uJn) =< 
Sp^ i~q, ~uJn)Sp- {q,ujn) >= uj'i+v%)qi whcrc v'^{q) = 
2pEVE{q) is the spin wave velocity. 

From the pole of the dynamic density-density correla- 
tion function, we can identify the speed of sound wave 
which is exactly the same as the spin wave velocity. This 
should not be too surprising. As shown in liquid '^He, 
the speed of sound is exactly the same as the phonon 
velocity. Here, in the context of excitonic superfluid, 
we explicitly prove that the sound speed is indeed the 
same as the spin wave velocity. From the analytical con- 
tinuation itUn Lo + we can identify the dynamic 
structure factor: S-{q,uj) = S-{q)S{uj — v{q)q) where 
S-{q) = pEqT^lv{q) is the equal time pseudo-spin corre- 
lation function shown in Fig.l. As g ^ Q^S-{q) q. 
The Feymann relation in BLQH which relates the dis- 
persion relation to the equal time structure factor is 



uj{q) = 



PETTq 



(6) 



which takes exactly the same form as the Feymann 
relation in superfluid '^He. Obviously, the Ve(<z) in 
the Fig. lb leads to the magneto-roton dispersion — 
q^VEiq) in the Fig.la. 
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Because the instability is happening at g = go instead 
of at g = 0, the vortex excitations in 0_ remain uncriti- 
cal through the transition. Integrating them out will not 
generate any singularities except the interactions among 
the pseudo-spin density 6p-. Expanding ¥5(9) near the 
minimum qo in the Fig. 1 leads to the quantum Ginsburg- 
Landau action to describe the ESF to the PSDW transi- 
tion: 



]5p- 



+ u{5p-f + v{5p-f + 



(7) 



where the momentum and frequency conservation in the 
Quartic and sixth order terms is assumed, A„ ~ ^ -i 
which is non-critical across the transition. While the 
corresponding quantity Ag ^ S- (q) in the phase rep- 
resentation EqnOis divergent, so Eqn [3] breaks down as 
q qQ and may not be used to describe the ESF to the 
PSDW transition. 

In sharp contrast to the conventional classical normal 
liquid (NL) to normal solid (NS) transition fv?\, the pos- 
sible cubic interaction term (Sp-^)^ is forbidden by the Z2 
exchange symmetry between the two layers dp- ^Sp^ . 
Note that the term in the first term stands for the 
quantum fluctuations of (5p_ which is absent in the clas- 
sical NL to NS transition. In the following section, we 
will show that the most favorable lattice is a square lat- 
tice. 

Lattice structure of the PSDW phase. In Eqnl?! 

r which is the gap of Ve((?) at the minimum tunes the 
transition from the ESF to the PSDW. In the ESF, r > 
and < (5p_ >= is uniform. In the PSDW, r < and 



< Sp- >=EG'^(^)e 



iG-x 



i(0) = takes a lattice struc- 



ture with reciprocal lattice vectors G. It was shown that 
in the classical NL to NS transition, due to the cubic 
term, a triangular lattice is the favorite lattice. However, 
due to the absence of the cubic term in Eqn[71 in the fol- 
lowing, we will show that the favorite lattice is the square 
lattice. At mean field level, we can ignore the uj depen- 
dence in Eqnlll Substituting < Sp^ Y.c^i^^^^'^^ 
into EqnlT] leads to: 
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FIG. 2: (a) The charge distribution of the PSDW in a square 
lattice. The "up " pseudo-spins take sublattice A, while the 
"down" pseudo-spins take sublattice B. (b) In general, there 
are quantum and thermal fluctuation generated vacancies in 
both layers denoted by a Q- 



Following [l7|, one can simplify Eqn[8]to: 



(9) 



where for the quadratic and quartic terms, all the con- 
tributions come from the paired reciprocal lattice vec- 
tors. After the scaling uq — > mT^I'^nQ, the quadratic 
term become the same for both lattices, but the quar- 



"G' 

he qr 

tic terms are still different Ua = 3(1 — l/m)M[l7|, so 



— 2^u,u/\ = 2^u. For the sixth order term, in 
square lattice, all the contributions are still from the 
paired reciprocal lattice vectors, however, in triangu- 
lar lattice, there is an additional contribution from a 
triangle where Gi + G2 + G3 = 0. After very care- 
ful counting, we find that = 5(3 — 9/m -I- 8/m'^)v 
from the paired contribution and v'^" = 5/6v from the 
additional contribution from the triangle in a triangu- 
lar lattice, so = 6jV,v/\ = 9|u. The mean field 
phase diagram of EqnlQ] is well known: If m > ( 
u < ), there is 2nd ( 1st ) order transition, there 
is a tri-critical point at u = 0. Minimizing / with re- 
spect to leads to — {—2ua + \/4u^ — 6var)/6va 
which holds for both Ua > and Ua < 0. Obviously, 
TT'O 7^ ^A- From Eqnini we can see that if u > 0, because 
ua < UA,va < VA, for any given n, fa{_n) < /a(»^), 
then /□(nn) < /□(?^A) < /a(»^a)- If m < 0, it is easy 
to show that when n'^ > nl = /□("-) < fA{n)- 

Because in the two lattice states, r < u\/2va, rig > 



2\u. 



2vu 



> 



MA 



9m 

' 6Sv 



> nt 



2+G3+(34,o this regime, we stiU have /□(nn) < /□(^a) < fAin^)- 
So we conclude that in any case, the square lattice is 



vyriA riA riA ua tia 5^ 

Z / Gi G2 G3 G4 Gs Ge Gi 



+ G2+G-i+Gi+G5+G6fi 



ttei favorite lattice. It has three elastic constants in- 



where u could be positive or negative, w > is always 
positive to keep the system stable. 

From EqnlSl we can compare the ground state energy 
of the two most commonly seen lattices: square lattice 
and triangular lattice. Square ( triangular ) lattice has 
TO = 4 ( m = 6 ) shortest reciprocal lattice vectors G = 
2'K I a where a is the lattice constant in a given layer. 



stead of two. Neglecting zero-point quantum fiuctua- 
tions, < bp- >= Ej 8{x — Ri) — J2i ^{^ — Ri — I) where 
the I = ^{ai+ 0,2) is the shift of the square lattice in the 
bottom layer relative to that in the top layer ( Fig. 2). 

This PSDW state not only breaks the translational 
symmetry, but also the Z2 exchange symmetry. It is 
very rare to get a 2d square lattice, because it is not a 
close packed lattice. Due to the special Z2 symmetry 
of BLQH, we show it indeed can be realized in BLQH. 
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The system is compressible with gapless phonon excita- 
tions determined by the 3 elastic constants. It is very 
interesting to see if in-plane soft X-ray or light scatter- 
ing experiments 0] can directly test the existence of the 
PSDW when dd < d < dc2- Note that the hght scatter- 
ing intensity may be diminished by a Debye- Waller factor 
due to the large zero-point quantum fluctuations in the 
PSDW 0. 

Vacancies, disorders and Coulomb Drag in the 
PSDW state. In principle, the Sp+ mode in Eqnl2] 
should also be included. It stands for the translational 
(or sliding ) motion of the PSDW lattice. However, any 
weak disorders will pin this PSDW state. Therefore, the 
(5/9+ mode can be neglected. Disorders will smear the 1st 
order transition from the ESF to the PSDW into a 2nd 
order transition. It was argued in [1] that disorders in 
real samples are so strong that they may even have de- 
stroyed the ESF state, so they may also transfer the long 
range lattice orders of the PSDW into short range ones. 
This fact makes the observation of the lattice structure 
by light scattering difhcult. The two square lattices in 
the top and bottom layer are locked together, so it will 
show huge Coulomb drag. However, vacancies generated 
by the large zero-point quantum fluctuations may play 
important roles in the drag ( Fig. 2b ). As the distance 
increases to the critical distance dd in Fig.l, the ESF 
turns into the PSDW whose lattice constant a — ^/inl is 
completely fixed by the filling factor i^i = 1/2, so it may 
not completely match the instability point l/qo- Due to 
this slight mismatch, the resulting PSDW is very likely 
to be an in- commensurate solid where the total number 
of sites Ns may not be the same as the total number of 
sites N even at T = 0. As the distance increases further 
del < d < dc2 in Fig.l, the PSDW lattice constant is still 
locked at a = ^/AttI. Assuming zero-point quantum fluc- 
tuations favor vacancies over interstitials [16[, so there 
are vacancies no even at T = in each layer. At finite 
temperature, there are also thermally generated vacan- 
cies na{T) ~ e~^^l'^ where A„ is the vacancy excitation 
energy. So the total number of vacancies at any T is 
ny{T) = no -\- na{T). Obviously, the vacancies in top 
layer are strongly correlated with those in the bottom 
layer. As shown in the correlated character of hop- 
ping transports of the vacancies between the active and 
passive layers can lead to a very large drag. We can es- 
timate the resistance in the active layer as i? ~ 1 /n^ (T) . 
As shown in [l8| , the drag resistance in the passive layer 
is Rd aoR where ao should not be too small. This 
temperature dependence is indeed consistent with that 
found in the experiment [Bj: starting from 200mK, Rd 
increases exponentially until to SOmif and then satu- 
rates. This behavior is marked different than that at 
both small and large distance where the system is in the 
ESF and FL regimes respectively. These vacancies also 
lead to the finite Hall drag in the PSDW. We conclude 
that in the presence of disorders, all the properties of the 



PSDW are consistent with the experimental observations 
in on the intermediate phase. The effects of very small 
imbalance on the ESN phase was investigated in ll| and 
was also found to be consistent with the experimental 
data in [l^ . 

Discussions. We compare the results achieved from 
the QGL theory with those achieved from the microscopic 
LLL-|~HF approximation in 6]. Especially, Cote, Brey 
and Macdonald in [6] numerically found that the lowest 
energy lattice structure of the PSDW is a square lat- 
tice. But it is not known if the HF+LLL calculations 
can describe the transition from the ESF to the PSDW. 
The QGL theory presented in this paper circumvents the 
difficulty associated with the unknown wavefunction at 
any finite d [l9| and treat both the interlayer and the 
intralayer correlations on the same footing, so can be 
used to capture competing orders on microscopic length 
scales and naturally leads to the PSDW as the interme- 
diate state which breaks translational symmetry. The 
theory puts the ESF state and the PSDW state on the 
same footing, characterize the symmetry breaking pat- 
terns in the two states by corresponding order parame- 
ters and describe the universality class of the quantum 
phase transition between the two states. We explicitly 
showed that the square lattice is the most favorite lat- 
tice. There are quantum fluctuation generated vacancies 
in the PSDW which lead to the unusual temperature de- 
pendence of the Coulomb drag observed in the exper- 
iment [Bj. The QGL theory is complementary to and 
goes well beyond the previous microscopic calculations. 
Recent experiments !20] found the real spin also plays 
important roles in some BLQH samples with large tun- 
neling gaps. Putting the spin into the theory remains an 
important open problem. 

J. Ye thank B. Halperin for helpful discussions and J. 
K. Jain for pointing out Ref. [l^ to us. 
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